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Different ways to incorporate two-dimensional systems, which are not 

amenable to separation of variables, into the framework of Supersymmetri- 

cal Quantum Mechanics (SUSY QM) are analyzed. In particular, the direct 

CN ■ generalization of one-dimensional Witten's SUSY QM is based on the super- 

s- 
charges of first order in momenta and allows to connect the eigenvalues and 



(N 



eigenfunctions of two scalar and one matrix Schrodinger operators. The use 
of second order supercharges leads to polynomial supersymmetry and relates 
a pair of scalar Hamiltonians, giving a set of such partner systems with al- 
most coinciding spectra. This class of systems can be studied by means of 
' new method of SUSY— separation of variables, where supercharges allow 

. separation of variables, but Hamiltonians do not. The method of shape 

| invariance is generalized to two-dimensional models to construct purely al- 

gebraically a chain of eigenstates and eigenvalues for generalized Morse po- 
Q_i' tential models in two dimensions. 



1. Introduction 



X. 

Supersymmetric Quantum Mechanics [1], [2] is an interesting framework to analyze non- 

relativistic quantal problems. In particular, it allows to investigate the spectral properties of a 

wide class of quantum models as well as to generate new systems with given spectra. SUSY QM 

gives new insight into the problem of spectral equivalence of Hamiltonians, which historically 

was constructed as Factorization Method in Quantum Mechanics [3] and as Darboux-Crum 

transformations in Mathematical Physics [4]. 

During last two decades SUSY QM became an important and popular tool to study a wide 

variety of quantum systems (see the list of reports presented to this Conference). It is easy to 

note that the main stream of the development in SUSY QM concerned one-dimensional models. 

Though the variety of multi-dimensional (especially two- and three-dimensional) problems is 

much wider and practically important, much less attention has been given in the literature to 
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the study of these models in SUSY QM. Thus the future progress seems to be mostly connected 
with investigation of its mult i- dimensional generalizations. The main aim of this paper (based 
on the talk at the Conference) is to summarize different results of previous investigations of 
two-dimensional SUSY QM. 

The paper is organized as follows. Two-dimensional generalization of the conventional 
Witten's formulation of SUSY QM is formulated in Section 2. In Section 3 the two-dimensional 
SUSY QM with the supercharges of second order in derivatives is presented. Section 4 contains 
a new, supersymmetric, approach for investigation of two-dimensional models, which are not 
amenable to separation of variables. This method is based on the second order supercharges 
introduced in the previous Section, and it gives a new opportunity to reduce the problem to the 
one-dimensional ones. Thus we obtain the specific method of SUSY— separation of variables. 
Section 5 is devoted to the generalization of the well known notion of shape invariance onto the 
two-dimensional models. Some new aspects typical for two-dimensional shape invariance are 
investigated. In Section 6 the deformation of SUSY QM algebra for models of Sections 3-5 is 
described. 

2. Direct two-dimensional generalization of the conven- 
tional Witten's SUSY QM 

The conventional one-dimensional SUSY QM was proposed by E.Witten [1]. It is charac- 
terized by the simplest realization of SUSY algebra: 

{Q + ,Q~} = H; (Q + ) 2 = (Q~f = 0; [H,Q k ] = 0, (1) 

where the superHamiltonian H is a diagonal matrix H = diag(h^°\h^), with = 
— <9 2 + VW(x), d = d/dx, and supercharges Q ± are off diagonal with elements q T = 
±d + dW(x) of first order in derivatives with superpotential W(x). In terms of components the 
(anti) commutation relations of SUSY algebra (1) mean, respectively, the factorization of Hamil- 
tonians, nilpotent structure of supercharges and intertwining of by q ± . The superpotential 
W(x) is defined by an arbitrary (possibly non-normalizable) solution *ff(x) = exp (—W(x)) of 
the Schrodinger equation h^fy(x) = e\&(x) with e < . If this solution is nodeless one has 
almost coinciding spectra of and or, equivalently, double degeneracy of the energy 
spectrum of H. 

The direct multi-dimensional generalization of the construction above was built in [5] both 
by direct extension [6] of one-dimensional formulas and using [7] the superfield approach of 
Quantum Field Theory. We call it as " direct" since it retains both all relations of superalgebra 
(1) and the first order form of the components of supercharges. 



Here we restrict ourselves to the particular case of two space dimensions x = (xi,x 2 ). 
The precise formulas for the 4x4 superHamiltonian and supercharges are then the following: 
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where two scalar Schrodinger operators , and 2x2 matrix Schrodinger operator are 
expressed in quasifactorized form in terms of the components of supercharges: 

= qfqf = -df + V^(x) = -df + (d^ix)) 2 - dfW(x); df = df + d 2 2 , 
= pjpl = -df + V {2 \x) = -df + (dtW(x)) 2 + dfW(x); 
h { S = Q^t + Pjvt = -hkdf + 5 ik WW(x )) 2 - dfW(x)) + 2d i d k W(x). 

The components of supercharges qf,pf are again of first order and depend on the two- 
dimensional superpotential W(x) : 

q± = T d l + d l W{x)- pf = e lk qt. (3) 

Two-dimensional 4x4 superHamiltonian and supercharges (2) realize the same conventional 
SUSY QM algebra (1). In components, the commutation relations in (1) are expressed as the 
intertwining relations between matrix and by operators q^,pf: 

h<£ q t = qrh<®; h$ Pk r=p-hM; = h^p 



+ 



(4) 



The energy spectra of and are in general different, but the intertwining relations (4) 
provide the equivalence of energy spectra between a pair of two scalar Hamiltonians h^°\ 
and 2x2 matrix Hamiltonian hff} . The "equivalence" means coincidence of spectra up to zero 
modes of operators qf,pf. Thus the supersymmetry (supersymmetric transformation) allows 
to reduce the solution of matrix Schrodinger problem with the Hamiltonian to solution of 
a couple of scalar Schrodinger problems h^°\ h^ 2 \ Due to the same intertwining relations the 
vector wave functions of matrix Hamiltonian h^) are also connected (up to a normalization 
factor) with the scalar wave functions of scalar Hamiltonians h^°\ : 

^ (1) (f;£) = q-¥°\x;E); < = 1, 2 ¥°\x; E) = q^ (x; E) 



W(x;E) = pr¥ 2 \x;E) ¥ 2) (x; E) = pf^ 1 ' (x; E). 

The Schrodinger operators with matrix potential are not something very exotic in Quantum 
Mechanics. In particular, the described two-dimensional generalization of SUSY QM was suc- 
cessfully used [8] to investigate the spectra of Pauli operator for fermion in external electromag- 
netic fields. Nevertheless, considering this rather good-looking construction, one question seems 



to be natural: is it possible to perform supersymmetric transformations in two-dimensional case 
avoiding any matrix Hamiltonians? 

3. Second order supercharges in two-dimensional SUSY 
QM 

The main idea, which could allow us to get rid of matrix components of superHamiltonian, is 
to explore the supercharges of second order in derivatives. For the first time such supercharges 
of higher orders in momenta were proposed for the one- dimensional situation in [9] (see also 
[10], [11]) leading to the polynomial deformation of SUSY algebra (see Section 6). In general, 
this approach implies a deformation of one relation of SUSY algebra (1) only, namely of the 
(quasi) factorization, but keeping unchanged the nilpotency of Q ± and the intertwining relations. 
The last seem to be the most important ingredient of SUSY methods in QM. 

The simplest variant of second order supercharges - of the so called reducible [10] form - 
gives uninteresting result in two-dimensional case: the intertwined partner Hamiltonians differ 
by a trivial constant only, and both of them admit the separation of variables (see details in [12], 
[13]). By this reason here we will be interested in general irreducible second order components 
of supercharges: 

q + = g lk (x)d t dk + Ciffldi + B(x); q~ = (q + )l (5) 

The important question we have to investigate now concerns the existence of Hamiltonians 

h (i) = _ A (2) + V {i) (x); i = 1,2; A (2) = dA, (6) 

which satisfy the intertwining relations with q ± of the form (5): 

h (D q + = q+h (2). q - h W = h <?) q - m (7) 

The first consequence of (7) restricts essentially the possible "metrics" gik(x) by dig ik + 
digik + d k gu = with solutions: 

#n = ax\ + aix-2 + h] g 2 2 = olx\ + a 2 xi + b 2 ; gn = -\ > {2ax 1 x 2 + a x xx + a 2 x 2 ) + b 3 . (8) 

This has to be taken into account in rewriting [12] the intertwining relations (7) in components: 

diC k (x) + d k Q(x) + A (2) g lk (x) ~ (V (1 \x) - V (2 \x))g tk (x) = 0; (9) 
A^Q(x) + 2d i B{x) + 2g lk (x)d k V^(x) - (V^(x) - V {2 \x))C t {x) = 0; (10) 
A^B(x) + g ik (x)d k d t V (2 \x) + Q(x)d t V {2 \x) - (V {1 \x) - V {2 \x))B(x) = 0. (11) 



The nonlinear system of second order differential equations (9) - (11) for unknown functions 
Ci(x), B(x), V^(x), V^(x) with constant parameters a,ai,bi in gik(x) does not admit the 
general solution, but one can look for its particular solutions with concrete metrics g^ and 
some ansatzes for unknown functions Cj. 

In particular, the system (9) - (11) is essentially simplified for the metrics of elliptic form 
gik{x) = 5ik- In this case all unknown functions in (9) - (11) can be found analytically in the 
general form [12], but for all such solutions both Hamiltonians h^ l \ turn out to admit 
the R— separation [14] of variables in parabolic, elliptic or polar c coordinates, i.e. this class of 
two-dimensional problems can be reduced to two one-dimensional models. 

Much more interesting situation appears for more complicated forms of metrics. Thus the 
list of particular solutions of the system (9) - (11) can be constructed analytically [12], [13], [15] 
for the hyperbolic metrics g^ = diag(+l, —1). Indeed, for this metrics a part of task of solution 
of (9) - (11) can be made in general form since it is reduced to a simpler system, which will be 
written in terms of coordinates x± = x\ ±x 2 d± = d/dx± and C + = C\ — C 2 ; C_ = C\ + C 2 . 
The general solution can be provided by solving the system: 

d.{C.F) = -d + (C + F); (12) 
d\F = d 2 _F, (13) 

where C± depend only on x±, respectively: C± = C±(x±). The function F, solution of (13), 
is represented as F = Fi(x + + x_) + F2(x + — xJ). The potentials V^'^{x) and the function 
B(x) are expressed in terms of Fi(2xi), F 2 {2x 2 ) and C±(x±), solutions of (12): 



V m2) = ±\(C + + C'_) + ^Cl + Cl)+ 1 -( K F 2 (x + -x^-F 1 (x + + x_)y, (14) 
\ (C+C. + F,(x + + x.) + F 2 (x + - , 
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where C means derivative in its argument. 

For lack of a regular procedure for solution of both equations of the system (12), (13), its 
particular solutions can be found starting from certain anzatses for functions C±(x±), F(x). 

1) Let C_ = 0, then from (12) one obtains F = <f)(x-)/C + (x + ). After inserting into Eq.(13) 
the separation of variables is possible, and particular solution reads d : 

^+(^+) = 7— 1— j 

5i exp(v A • x+) + 5 2 exp(— v A • x + ) 

F\ )2 (2x) = 5i(7i j2 exp(2v / Ax) + 5 2 <y 2yl exp(-2v / Ax), 

the Greek letters - arbitrary constants, depending on sign of A they may be real/complex. 

c Thc reducible second order supercharges just correspond to separation of variables in polar coordinates. 
d Here and below are defined only up to an arbitrary real constant: F\ — > F\ + e, F 2 — > F 2 — e. 



2) Let F(x) allows also the factorization: F = F + (x + ) ■ F_(x_). Then from Eq.(12): 

x± 

C ± = ^±^-J F ± dx' ±: (15) 

and there are two options to fulfill the condition (13), i.e. F(x) = Fi(2x\) + F 2 (2x 2 ): 

a) F±(x±) — e±x±, b) F± = a± exp(\/\ • x±) + 5± exp(— v^A • x±). (16) 

Corresponding potentials can be found according to Eq.(14), being similar to ones obtained 
in [16] in quite different approach. Below some other solutions of (12), (13) will be built [15]. 

3) Let us start now from the general solution of (12): 

dx-i. f dx- 1 



-)/(C + C_). (17) 



C + J C_ 

Then Eq.(13) gives the functional-differential equation for the functional L(A + — A J) with 
A' ± = 1/C ± (x ± ) : 
/A'" A"'\ 

- jr)L(A + - A_) + 3« + A"_)L\A + - A_) + {A* - A'1)L"{A + - A_) = 0, (18) 

where L' denotes the derivative of L with respect to its argument. If we take functions A± 
such that A" ± = \ 2 A±, A = const, Eq.(18) will become ordinary differential equation for L with 
independent variable (A + — AJ). It can be easily solved: 

L{A + -A_) =a{A + - A_y 2 + (3, 

where A± = <j±exp{\x±) + 5±exp(— \x±) with a + • 5 + = o"_ ■ 5_ and a, (5 - real constants. For 
A 2 > 0, choosing a± = —5± = k/2 or a± = +5± = k/2, we obtain (up to an arbitrary shift in 
x±) two particular solutions: 

3a) A± = ksmh(\x±), 36) A± = k cosh(\x±). 

Then (17) leads to: 

3a) F 1 (2x) = — -fy— + k 2 cosh(2Ax); (19) 
cosh (Ax) 

36) F l( 2,) = -F 2 (2,) = -^ + * 2 rf(A,), C ± = --±- y k^O. (20) 

For A 2 < hyperbolic functions must be substituted by trigonometric ones. 

We have to remark that the case A 2 = 0, i.e. A'± = 0, is not of interest, leading to 
trivial superpartners. However, choosing in (20) A — > 0, k, k±, k^ 1 — > simultaneously, so that 
A 2 ~ &i ~ k 2 l ~ A; 2 , we obtain the solution: 

~ k 
Fi(2x) = -F 2 (2x) = hx~ 2 + k 2 x 2 , C± = — . (21) 

x± 



One can check that (12) is also satisfied by 

F 1 (2x) = -F 2 (2x) = k lX 2 + k 2 x A , C± = ±—. (22) 

x± 

4) Starting again from (17), it is convenient to pass on to new variable functions C± = 
±/±//±. Then F in (17) is represented in the form F = U(f + f)f' + f'_ with an arbitrary 6 
function U. After substitution in (13) one obtains the functional-differential equation: 

/ f" f" \ / f" f" \ 

(tff- - flf a )u"U) + 3/(^ - J T )u'U) + {jr- jr)U(f) = o, / = /+/_. 

For particular form of functions f± = a±exp(Xx±) + fl±exp(—Xx±), this equation becomes an 
ordinary differential equation for U with independent variable /. Its solution is U — a + 4bf + f- 
(a, b— real constants). Then functions 

F 1 (x) = ki(a + a_ exp(Ax) + (3 + {3- exp(— Xx)) + k 2 (a 2 + a 2 _ exp(2Ax) + P+fi 2 . exp(— 2Ax)), 

— F 2 (x) = ki(a + P- exp(Ax) + (3 + a- exp(— Ax)) + k^a+fi 2 , exp(2Ax) + f3\a 2 _ exp(— 2 Ax)), 
± a± exp(Ao: ± ) + (3± exp(-Ax ± ) 
X(a±exp(Xx±) — P±ex.p(—Xx±)) 

(with ki = aX 2 , k 2 = 46A 2 ) are real solutions of (12), (13), if a±, (3± are real for A 2 > 0, and 
a± = (3* ± for A 2 < 0. 

5) To find a next class of solutions it is useful to rewrite (12) in terms of x± i2 : 

(Fx(2a;i) + F 2 (2x 2 ))d 1 (C + + C_) + F[(2 Xl )(C + + C_) + F! 1 {2x 2 ){C + - C_) = 0. 
Among known particular solutions the most compact one is: 



6 



2 



C+(a;) = C_(x) = ax 2 + c, Fi(2xi) = 0, F 2 (2x 2 ) = —. (24) 

x 2 

After inserting these solutions (19) - (24) into the general formulas (14), one obtains the 
analytical expressions for potentials. Their explicit form can be found in [15]. 

The additional class of particular solutions of the system (9) - (11) obtained for the case of 
degenerate metrics = diag (1,0) can be found also in [15]. 



4. SUSY— separation of variables 

From the very beginning in this paper we are interested in two-dimensional quantum systems, 
which are not amenable to separation of variables. The supersymmetric approach, namely the 
intertwining relations (7), allows to formulate some specific supersymmetric alternative to the 

e Due to Eq.(13), the function F should be additionally representable in the form F = Fi(2xi) + ^2(2x2). 



conventional notion of separation (including the so-called R— separation [14]) of variables. The 
main idea [17] is to consider such particular class of solutions of intertwining relations (7), when 
the components of supercharge q* 1 are amenable to separation of variables but Hamiltonians 
are not. In this case the Hamiltonians turn out to be partially solvable, or in another 
terminology, quasi-exactly-solvable [18]. Both terms above mean that a part of spectrum (and 
possibly of corresponding eigenfunctions) of Hamiltonian is known. The crucial ingredient of 
the approach [17] is in the investigation of zero modes of intertwining operators q ± . 

The general scheme of the method is the following. Let us suppose that N + 1 normalizable 
zero modes of q + are known (for example, due to separation of variables in q + ): 

q + a n (x)=0; n = 0,1, ...,7V; q + (l(x) = 0, (25) 

where Q(x) is a column vector with components Q n (x). From the intertwining relations (7) one 
can see that the space of zero modes is closed under the action of hP\ and therefore: 

h (2) Q{x) = CQ{x), (26) 

where (7 = | |cj fc | | is a c— number x— independent real matrix. If the matrix C can be diagonalized 
by a real similarity transformation: 

BC(B)- 1 = A = diag(\ , X u X N ), (27) 

the problem is reduced to a standard algebraic task within the zero modes space: 

h (2 \M{x)) = A(M(x)). (28) 

It is not clear in advance, whether this general scheme is realized practically? To put it 
differently, are there any solutions Cj(x), B(x) of the intertwining relations (9) - (11), which 
give q + with separation of variables? 

To investigate this problem, it is useful [17] to transform the supercharge q + by the special 
similarity transformation, which removes the terms linear in derivatives: 

q+ = e (-x(*)) 5 + e (+x(50) = 9 2_ 9 2 + l( j p 1 (2x 1 )+F 2 (2x2)); x{x) = ~( J C+{x + )alx + + J C-(x-)dx- 

(29) 

These new operators q + obviously obey the condition of separation of variables realizing the 
first step of our scheme of SUSY— separation of variables. Zero modes of q + can be found as 
linear superpositions of products of one dimensional wave functions r) n (xi) and p n (x 2 ), satisfying 
Schrodinger equations (with e n - the separation constants.): 

(-<9i - ^-Pi(2a;i))77 n (a;i) = e n r/ n (a; 1 ); (-d 2 + ^F 2 {2x 2 )) pn{x 2 ) = e n p n {x 2 ). (30) 



In analogy to (29), one can define operators 

h = exp(- x (f)) exp(+ x (f)) = -df + C^x)^ - C 2 (x)d 2 - ^F 1 (2x 1 ) + ^F 2 (2x 2 ), (31) 
and eigenfunctions of q + as: 

tt n (x) = exp(-x(£)) • fi„(x), (32) 

keeping however in mind that the normalizability and orthogonality are not preserved auto- 
matically due to non-unitarity of the similarity transformation. 
Then using (30) one can write: 

m n (x) = [2e n + Ci(£)3i - C 2 (x)d 2 ]Q n (x). (33) 

It is not evident from (33), but the space spanned by functions fl n {x) is closed under the action 
of h. It will be demonstrated explicitly in the concrete model below. 

In contrast to (29), where variables are separated, no separation for h, which would make the 
two-dimensional dynamics not-trivially reducible to one-dimensional dynamics. In this regard 
we refer to this method [17] for partial solvability as to SU SY— separation of variables. 

The scheme of SUSY— separation of variables formulated above can be used for arbitrary 
models satisfying the intertwining relations (7). The list of such models is already rather long, 
and it may increase in future, but it is very important to check the applicability of the scheme 
on the concrete model where the explicit solutions can be constructed. Actually, it means that 
solutions of two one-dimensional problems (30) can be found analytically. Below we briefly 
describe such a model - generalized two-dimensional Morse potential. 

Among the solutions [15] of the system (12), (13) we focus attention on the particular case 
with a specific choice of parameters and A > 0, a > 0, a-real constants ^ : 

(XX 

C+ = 4aa; C_ = 4aa ■ coth — (34) 
fi(xi) = -J l {2x l ) = -A(e- 2 ^-2e-^- 2=1,2; (35) 



V m2) = a 2 a(2a t 1) sinh" 2 (^-) + 4aV + A 



^-2axi _ 2e~ axi -\- e ^ 2ax 2 _ 2e~ aX2 



(36) 



One easily recognizes in (36) a sum of two Morse potentials plus a hyperbolic singular term 
which prevents to apply the conventional methods of separation of variables. These singular 
terms can be both attractive, for \a\ > |, or one repulsive and one attractive, for \a\ < \. The 
parameter a will be further constrained by the condition that the strength of the attractive 
singularity at x^ — > should not exceed the well known bound — l/(Ax 2 __ ). 
•^For the complexification of the model see [19] 



The normalizable functions Cl n (x) (and fl n (x)) can be constructed from the well known [20] 
normalizable solutions of (30) with e n < : 

Q n (x) = exp(-^^)(66) s "i 71 (-n, 2s n + 1; &)F(-n, 2s n + 1; 6), (37) 

where F(—n, 2s n + 1; £) is the standard degenerate (confluent) hypergeometric function, reduc- 
ing to a polynomial for integer n, and 



2^ 1 

exp(-axi); s n = n - - > 0; e„ = -A 

a a: z 



a , 1 , 



1 2 



1 "71 (T ' + 2 ) 



(38) 



The number (N + 1) of normalizable zero modes (37) is determined by the inequality s n > 0. 

The condition of normalizability of zero modes fl n (x), together with the absence of the "fall 
to the centre" , leads [17] to the following two ranges of parameters: 



1 1 v A 1 
a E (— oo, ;=); s n = n > —2a > 0. (39) 

oG(-ii); s >2(|a| + l). (40) 

In Sections 4 and 5 only the region (39) will be considered. Inequalities (39) can be satisfied 
by the choice of a and A, and/or by suitable restriction on the number N of zero modes fl n (x). 

Analysis of the action of h in (33) gives that the matrix C is of triagonal form. It can be 
diagonalized explicitly by a similarity transformation, and the eigenvalues of coincide 
with its (all different and nonzero) diagonal elements: 

E k = Cfcfc = -2(2aa 2 s k - t k ). (41) 

The resulting eigenfunctions of are obtained (see Eqs.(27), (28)) from the constructed zero 
modes fl n {x) and the similarity transformation B : 

^! N - n {x) = EfLobMffl. (42) 

For the algorithm of iterative construction of coefficients b n i see [17]. Thus the construction of 
the set of eigenfunctions, which lie in the space of zero modes, is completed. 

These eigenfunctions ^(af) may be also used for constructing more general eigenfunctions 
of via a product ansatz: 

$(x) = tt fc (x).6(x). (43) 

Three such eigenfunctions based on \l>o were constructed in [17]. Within the bounds imposed 
(39) only one of them is normalizable, though for the region (40) all three are normalizable. 



5. Shape invariance in two dimensions 



In the previous Section we developed the method which led to construction of the partially 
solvable (quasi-exactly-solvable) two-dimensional models. Let us remind now the well known 
in one dimension and very elegant method of shape invariance [21] usually associated with 
the exactly solvable one-dimensional systems. Our aim here is to generalize the idea of shape 
invariance to the two-dimensional case [17]. 

Refering readers to the original paper [21] and reviews [2] for the detail discussion of standard 
one-dimensional shape invariance method, let us list its main steps only: 

H(x;a) = H(x;a) + 71(a); a = a(a) (44) 

where 7Z(a) is a (c-number) function of a. The absence of spontaneous breaking of supersym- 
metry for all values of a implies that the lowest eigenvalue Eq(o) of H(a) vanishes and the 
corresponding eigenfunctions ^o(a) are normalizable zero modes of Q + (a). 

The intertwining relations Q~(x; a)H(x; a) = H(x; a)Q~(x; a) with the standard first order 
supercharge allow in this case to solve the entire spectral problem for H(x; a). The crucial steps 
are as follows. 

H(x; a)^ (x; a) = E (a)^ (x; a) = 0; H(x; a)^ (x; a) = 7t(a)^ (x; a). (45) 

It is important to remark that ^ (x;a) = ^ (x;a) has no nodes and therefore is the ground 
state wave function of H(x; a). Then 



Provided 



H(x; a) 



Q (xa)^/ (x; a) 



Q (x;a)^ (x;a) 



= 71(a) 



Q (x; a)ty (x; a) 



(46) 



is normalizable, we have generated an excited state of H(x; a), and 
thus 71(a) is naturally positive. It is clear that these steps can be repeated up to the last one, 
where the resulting wave function ^ will no more be normalizable. 

It is also clear that the isospectrality of H(x; a) and H(x; a) (up to the only zero mode 
^/q(x;o)) implies that there is no eigenvalue of H(x;a) between zero and the ground state 
energy E (a) of H. This observation leads to a proof that after suitable iterations one gets 
the entire spectrum of H(x;a). This method is referred as algebraic solvability (or complete 
solvability) by shape invariance in one-dimensional SUSY QM. 

To proceed to the formulation of two-dimensional shape invariance, we start from the rela- 
tively simple two-dimensional case of systems with conventional separation of variables: 

H(x) = H^xi) + H 2 (x 2 ); x=(x 1 ,x 2 ). 

Now suppose that Hi and H 2 both are shape invariant: 

Hi(xi] a,i) = Hi(xi, ai) + 7li(ai) <-> H(x ; a) = H(x ; a) + 71(a); a = (a ± , a 2 ). (47) 



In order to realize a nontrivial intertwining relations for H, H one can consider factorized 
supercharges of second order written as products of first order supercharges: 

Q ± = Qt-Qh Qf = Td t + W l (x l ). (48) 

There is the considerable difference with respect to the one-dimensional case. The crucial 
reason is that the space of zero modes of supercharges becomes now of higher dimensionality 
including the products of one- dimensional zero modes of the first Hamiltonian times all states of 
the second Hamiltonian and vice versa. While iterations are again obviously possible, it is clear 
that one can not argue about the entire solvability of the spectral problem, because in general 
many zero modes of (48) exist. Their number depends on the confining properties of Hi and H2. 
For example, in a case of oscillator-like potentials this number becomes infinite, and they are 
distributed over the whole spectrum. In this case only partial solvability of H can be achieved 
by the choice of (48) and shape invariance. Of course, one can solve such trivial models by 
separate use of Q^ 1 = Qf, which allows to solve the entire spectrum of two-dimensional model 
in terms of the one-dimensional ones. 

Let us suppose to have two-dimensional system (without separation of variables) with a 
Hamiltonian H, which is related to H by (44). For simplicity (in general, there is no connection 
between the dimensionality of the Schrodinger equation and the dimensionality of the parameter 
manifold), we assume that shape invariance is realized with one parameter a. Two-dimensional 
SUSY QM does not identify here zero modes of Q ± with the ground state of the Hamiltonian. 
Thus one has to repeat the steps (44) - (46) above by taking into account E (a) 7^ 0. In order 
to make our discussion more explicit we will from now on refer explicitly to the model (34) - 
(36) (H = h^), H = with the parameter a being bound to (39). 

First of all we observe that this model is indeed shape invariant (the infinite domain given 
by (39) allows iterations of (44)): 

a = a - \] 11(a) = a 2 {Aa - 1)) (49) 



The starting point is to write (46): 



H(x; a) 



(E (a-±)+K(a) 



Q (x;a)^ (x;a- i) 



(50) 



where E (a) and ^o(x; a) not to be identified with ground state. Thus we have constructed the 
new eigenstate and eigenvalue of H(x; a), provided Q~(x; a)^/ (x; a — |) is normalizable. Note 
that the eigenvalue (^E (a — \) + TZ(a)j is larger than E (a) with the bounds of (39). 

It is interesting that the energy of the first iteration of shape invariance in (50) coincides 
precisely with the eigenvalue of the additional solution <3>(x) mentioned in the previous Section: 

E = E {a - 1) + K(a) = a 2 [4a(l - s ) + (2s - 1)] + 2e . 



H(a) 



Q-(a)Q-(a-^ (a-l) 



Q-{a)Q-{a- l -)^! Q {a-l) 



The next iteration of shape invariance will give: 

£ (a-l)+ft(a-i)+ft(a)) 

(51) 

and the new eigenfunction Q~(a)Q~(a — |)\E'o( a — 1) can be written explicitly as function of x. 
Provided normalizability is ensured, one can thereby construct a chain by successive iterations 
of (50) and (51), since Q~{a) has no normalizable zero modes in (39). The end point of such a 
chain will be given by non-normalizability of the relevant wave function. 

Let us stress that though we illustrated both methods by the same model, SUSY— separation 
of variables can be implemented completely independently from shape invariance. For example, 
the model considered in the range (40) admits the method of SUSY— separation of variables, 
but shape invariance can not be applied since the domain (40) is too small. 



6. Polynomial algebra for two-dimensional SUSY QM 

In this last Section we will analyze the deformation of conventional SUSY QM algebra 
(1) due to introducing of second order components of supercharges in Sections 3-5. It is 
obvious that keeping the intertwining relations (7) and the matrix structure of supercharges 
Q ± we do not change two relations of SUSY algebra (1), which express the supersymmetry of 
H and nilpotency of Q ± . But the third relation, which gives the (quasi)factorization of the 
components of H, cannot be fulfilled. In the one-dimensional case the anticommutator of Q ± 
gives [10] the operator of fourth order in derivatives which can be represented as a second order 
polynomial of the superHamiltonian H. The situation changes crucially for two-dimensional 
systems of Sections 3-6, where in general a new diagonal operator of fourth order appears: 

R={Q + ,Q~} (52) 

This operator obviously commutes with the superHamiltonian H due to supersymmetry of H 
(intertwining relations (7)). It is shown in [13] that for Laplacian metrics in (5) g^ = 5n~ (where 
variables can be separated) operator R can be reduced to the second order symmetry operator R 
up to a second order polynomial (with constant coefficients) of H. But for all other metrics g^, 
including hyperbolic and degenerate ones, the Theorem was proved [13]: The symmetry operator 
R is essentially of fourth order in derivatives, its order can not be reduced. The components 
i = 1,2 of R are the symmetry operators of the components : [h^\ R^} = 0. The 
explicit expressions of R^ for known solutions of (9) - (11) can be written (see [22]). Thus all 
two-dimensional models of Section 3 are integrable 9 . 

s We note that integrability of the system does not mean in general its (even partial) solvability. 



As for generalized Morse model of Section 4, the quantum integral of motion R = Q~Q + 
gives zero on the eigenfunctions ^(x) by construction, since they are zero modes of Q + . But by 
a direct calculation one can check that all three additional eigenfunctions <&(x) (see (43)) of H 
are simultaneously [17] eigenfunctions (with nonzero eigenvalues) of the symmetry operator R. 
Thus they belong to a system of common eigenfunctions of two Hermitian mutually commuting 
operators H and R. 
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